ASIAN JOURNAL OF MATHEMATICS AND APPLICATIONS 
Volume 2014, Article ID ama0114, 11 pages 

ISSN 2307-7743 

http://scienceasia.asia 


SOLUTION OF THIRTEENTH ORDER BOUNDARY VALUE 
PROBLEMS BY DIFFERENTIAL TRANSFORMATION METHOD 


MUZAMMAL IFTIKHAR, HAMOOD UR REHMAN, AND MUHAMMAD YOUNIS 


ABSTRACT. The aim of this paper is to use the differential transformation 
method (DTM), a semi-numerical and semi-analytic technique for solving lin- 
ear and nonlinear thirteenth order boundary value problems. The approximate 
solution of the problem is calculated in terms of a rapidly convergent series. 
Two numerical examples have been considered to illustrate the efficiency and 
implementation of the method and the results are compared with the method 
developed in [1]. 


1. INTRODUCTION 


Higher order boundary value problems occure in the study of fluid dynamics, 
astrophysics, hydrodynamic, hydro magnetic stability, astronomy, beam and long 
wave theory, induction motors, engineering and applied physics. The boundary 
value problems of higher order have been examined due to their mathematical im- 
portance and applications in diversified applied sciences. 

The theory of differential transformation was established by Zhou [10]. The au- 
thor solved linear and non-linear IVPs arising in circuit analysis. Pukhov [4] also 
studied differential transformation method at the same time. The method con- 
structs an analytical solution for differential equations in the form of a series. It is 
a semi-numerical and semi-analytic technique that formulizes the Taylor series in 
a totally different manner. It differs from the traditional higher order Taylor se- 
ries method, which computationally, takes more time for higher orders BVPs. The 
DTM transforms the given BVP into a recurrence relation that finally leads to the 
solution of a system of algebraic equations as coefficients of a power series solution. 
The method is useful to obtain both exact and approximate solutions of linear and 
non-linear BVPs. There is no need to discretization, linearization or perturbation, 
large computational work and round-off errors are avoided. 

Siddiqi and Iftikhar used the variation of parameter method for solving the seventh- 
order boundary value problems in [7]. Siddiqi and Akram [5] used nonic spline 
and non-polynomial spline technique for the numerical solution of eighth-order lin- 
ear special case boundary value problems. These methods have also been proven 
to be second order convergent. Recently, Akram and Rehman presented the nu- 
merical solution of eighth-order boundary value problems using the reproducing 
kernel space method [2]. Geng and Li [3] construct a reproducing kernel space 
and solve a class of linear tenth-order boundary value problems using reproducing 
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kernel method. Siddiqi et al. [6] used the variational iteration technique for the so- 
lution of tenth order boundary value problem. Siddiqi and Akram [8] presented the 
numerical solutions of the tenth-order linear special case boundary value problems 
using eleventh degree spline. Siddiqi and Iftikhar [9] determined the approximate 
solutions of seventh, eighth, and tenth-order boundary value problems using homo- 
topy analysis method (HAM). Adeosun et al. [1] presented the variational iteration 
method (VIM) to find the approximate solutions of linear and nonlinear thirteenth 
order boundary value problems 

In the present paper, thirteenth-order boundary value problems are solved using 
DTM. The following thirteenth-order boundary value problems are considered 


u03) (x) = f(x, ,u(x)), a<xadd, 
(1. 1) u%(a) = Ai, 
uD(b) = Bj 


where for i = 0, 1, 2, ...,6 and j =0,1,...,5. Aj’s and B;’s are finite real constants. 
Also f(x, u(x)) is a continuous function on [a, b]. 


In the following section the differential transformation method is discussed. 


2. DIFFERENTIAL TRANSFORMATION METHOD 


The differential transformation of the kth derivative of a function f(x), at x = £o, 
is defined as 


(2. 1) F(k) = a | 





dak 


of) 


and the inverse differential transformation of F'(k) is defined by 


F(k)(a — zo)". 


M: 


(2. 2) f(x) 


In real applications, the function f(a) can be expressed as a finite series and Eq. ( 
2. 2 ) can be written as 


(2. 3) f(x) 


N 
So FO) = zo)". 
k= 


0 


Substituting F'(k) from Eq. ( 2. 1 ) into Eq. ( 2. 2 ), gives 





oS 1 [dE f(x 
(2. 4) fe) =e- a || 


dak 
k=0 


which is the Taylor series for f(x) at x = xo. From Eq. ( 2. 1 ) and Eq. ( 2. 2), 
the following theorems can be proved. 
Theorem 1 If f(x) = g(x) + h(x), then F(k) = G(k) + H(k). 
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Proof Using the definition ( 2. 1 ), yields 








(2. 5) ow = 5 (|. 

(2. 6) H(k) = a [| 

(2 7) Fi) = [gale tro 
Using Eqs. (2. 5) - (2. 7), yield 

(2. 8) F(k) = G(k) + H(h). 
Similarly, it can be proved that 

(2. 9) F(k) = G(k)— H(h). 


Theorem 2 If f(x) = cg(x), then F(k) = cG(k), where c is a constant. 
Proof Using the definition ( 2. 5 ), 








(2. 10) cw = [| 
Fk) = 4 [irten] 

(2. 11) = cn Tao 

Using Eqs. ( 2. 10 ) and ( 2. 11 ), give 

(2. 12) F(k) = cG(k). 


Theorem 3 If f(x aoe) then F(k) = (k+m)G(k+™m). 


) = 
Proof By definition ( 2. 5 ), 











© = alae S 
as 2 a (k+m) ee 
then 
(2. 14) F(k) = EEE ack +m) 
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Theorem 4 If f(x) = g(x)h(x), then F(k) = Yp, -o G(k1)H(k — kı). 
Proof By definition ( 2. 5 ), 


(2. 15) F(0) = 


lr Sie 





ll 

| es | 

a 

aJa 

Ble 
EN 
= 
Aisa 
A 
tea 

a 

S 

) EE j 

8 

i 


(2. 16) z 
FQ) = 5 geno] 


(2. 17) = G(2)H(0)+G(i)H D + GOH(2). 


Q 
E 
x 
© 
+ 
Q 
€ 
= 


In general, it can be written as 


k k2 


(2. 18) F(k) = X X G(k1)H (kz — ki) L(k — k2). 


k2=0 kı =0 


Theorem 5 If f(x) = eò, then F (k) = Ai 
Proof Using definition ( 2. 5 ), gives 


F(k) = kl Fz | a 
AF 


(2. 19) = 7 


RT. 
2 


Theorem 6 (i) If f(x) = Sin(az + b), then ay) 
os (2 + 


Esin (F +0). 
(ii) If f(x) = Cos(ax + 8), then F (k) = B). 


Proof (i) Using definition ( 2. 5 ), gives 


F(k) = L zl l Asins + 8) 
dx zeg 
k 
(2. 20) = Sin (= J B) : 
(ii) Using definition ( 2. 5 ), gives 
1 k 
F(k) = A È zCos(az +8) a 
(2. 21) = a cos (F + B) : 


Theorem 7 If f(x) = g(x)h(x)l(x), then 
F(k) = Eka- Dr-o (kx) H (ke — kı )L(k — ka). 
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Proof By definition ( 2. 5 ), 




















(2. 22) FO = A WEEE] = G(0)H(0)L(0) 
Fa) = 7 [loner] 
= Enea) g T ie + alan) | 
(2. 23) = G(1)H(0)L(0) + G(0)H(1)L(0) + G(0) H(0) L(A), 
Fe) = 5 RO] 
= G(2)H(0)L(0) + GA) H(1)L(0) + G1) H(0)L(1) 
+G(1)H(1)L(0) + G(O)H(2)4(0) + GOHO) 
(2. 24) 4G(1)H(0)L(1) + G(0)H(1)L(1) + G(0) H(0)£(2). 
In general, it can be written as 
k ko 
(2. 25) F(k) = X X G(k1)H (kz — ki) L(k — k2). 
k2=0 ki=0 


2.1. Analysis of the Method. The following nth order BVP is considered 
(2. 26) ul”) = f(x, u, uu ys uD), ax<x<bd, 


with the boundary conditions 








du d?u d°VDy 
2.2 B vee ; 
The differential transformation of the problem ( 2. 26 ) is 

F(k) 
22 U(k = 
(2. 28) (ktn) = Gy, 
where F(k) is the differential transformation of f(x, u, w, u”, ,u("-)). 
The transformed boundary conditions can be written as 
N m-1 

(2. 29) U(k) = A,U(m) = (k—i)U(k) = Bm, (m <n), 

k=0 i=1 


where m is the order of the derivative in the boundary conditions and A, Bm are 
real constants. 

Using Eqs. ( 2. 28 ) and ( 2. 29 ), the values of U(i), i = 1,2,3,--- can be 
determined and then using the inverse differential transformation, the following 
approximate solution upto O(a%*+!) can be determined 


N 
(2. 30) t= 5S s"U(k). 
k=0 


To illustrate the applicability and effectiveness of the method, two numerical ex- 
amples are considered in the following section. 
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3. NUMERICAL EXAMPLES 


Example 3.1 The linear thirteenth order BVP, is considered as 
(3. 1) u(x) = Cosa — Sina, 0<a<1, 


subject to the boundary conditions 


u(0) = il, u(1) = Cos(1) + Sin(1), 
vw)(0) = 1, u®(1) = Cos(1) — Sin(1), 
u?)(0) = -1, u (1) = —Sin(1) — Cos(1), 

(3. 2) u8(0) = -1, u®(1) =—Cos(1) + Sin(1), 
u4(0) = 1, uM (1) = Cos(1) + Sin(1), 
u)(0) = 1, wu (1) = Cos(1) — Sin(1). 
w(0) = -1, 


The exact solution of the problem is u(x) = Cosx + Sinz. 
Applying the theorems 1, 3 and 6 the differential transformation of the problem ( 
3. 1 ) can be determined as 





k! wk 1 wk 
. k+1 - ) i 
(3. 3) U(k +13) (k+ 13)! non (=) asin (=) 
Using the Eq. ( 2. 1 ), the boundary conditions ( 3. 2 ) at a = 0, can be 
transformed as 
ey UO) = LUN) =L) =, A 
U (3) = 5 U(A = 7, U (5) = 3, U(6) = =; 
Xio U(k) = Cos(1) + Sin(1), 


( 
Sono KU(k) = Cos(1) — Sin(1), 
(3.5) Zgo K(k = IU (k) = —Sin(1) — Cos(1), 





(k— 1)U(k) = 

Yeo klk — 1)(k — 2)U(k) = —Cos(1) + Sin(1), 
Si o klk — 1)(k — 2)(k — 3)U(k) = Cos(1) + Sin(1), 
Seg klk — 1)(k — 2)(k — 3)(k — 4)U(k) = Cos(1) — Sin(1), 


where n is a sufficiently large integer. Using the recurrence relation ( 3. 3 ), the 
transformed boundary conditions ( 3. 4 ) and the inverse differential transformation 
( 2. 3 ), the following series solution for n = 13, can be written 








r2 r3 xt gp? 
| A B Cx? 
u(x) i te ret ag aon am z+ Bx? + Cr 
x!3 
.6 De" + Ext! + Fr? 4 + O(x14 
3:8) ee eee w+ Sea7on0800 tO”): 


where, according to the definition ( 2. 1 ) 





ee ee mn ae ama 
(16) (11) u2) 
“e =U(10), E © = y(11), F = *5® =U(12). 
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From the transformed boundary conditions ( 3. 5 ), the following system of linear 
equations in terms of A, B, C, D, E and F can be determined 























8605396801 
6227000800 TAT BTC + D+ E+ F=Cos(l] + Sinfl], 

143700479 

Fogo + 74t 8B +9C + 10D +11E + 12F = Cos[l] — Sinfi], 

54885599 | 

3001Gs00 T 424 + 56B + 72C + 90D + 110E + 132F = —Cos|l] — Sin[1], 
1209601 910A + 336B + 504C1 + 720D1 + 990F + 1320G = -C [1] + Sinfi] 
3628800 ` i = —Cos 
544321 
360380 T 8404 + 16808 + 3024C1 + 5040 D1 + 7920F + 11880G = Cos|1] + Sin[1], 

1 
49399 + 25204 + 6720B + 15120C1 + 30240D1 + 55440F + 95040G = Cos|1] — Sin[1]. 


The solution of the above system, yields 

A = —0.00019841261645169112, B = 0.00002480111056133196, 

C = 2.7568717766803008 x 1078, D = —2.769904021536794 x 1077, 
E = —2.4115063057288555 x 1078, F = 1.8105803100768588 x 107°. 
The series solution can, thus, be written as 


r2 r? rt x5 xê 








= 1 ' ' 0.0001984132" + 0.00002480112° 
u(x) t-37 6 + 54+ o 720 ce £ 
+2.75687 x 107x? — 2.7699 x 1077 x1? — 2.41151 x 10782! + 1.81058 x 10792!” 
gi3 i 
(3.8) + §a570R0800 0C >) 


The comparison of the approximate solution of Example 3.1 obtained with the help 
of DTM and the approximate solution using VIM [1] is given in Table 1. From 
the numerical results, it is clear that the DTM is more efficient and accurate. By 
increasing the order of approximation more accuracy can be obtained. The results 
are also expressed graphically in Figure 1. The solid line represents the curve 
corresponding to the exact solution whereas the dotted line corresponds to the 
approximate solution. In Figure 2, absolute errors are plotted. 
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FIGURE 1. Comparison of 


the 


approximate 


solution 


with the exact solution for 


problem 3.1. 
approximate solution, 


Dotted line: 
solid 


line: the exact solution. 








FIGURE 2. Absolute errors 
for problem 3.1. 


Example 3.2 The thirteenth order non-linear BVP is considered as 


(3. 9) 


subject to the boundary conditions 


(3. 10) 


u(x) = eu? (x), O<a< 
u(0) = 1, u(1) =e, 
wo) = 1, wh) =e, 

u® (0) = 1, u@(1)=e, 
u®(0) = 1, u®(1)=e, 
u® (0) = 1, uM(1) =e, 
u®(0) = 1, «®()=e 

u® (0) = 1, 


The exact solution of the problem is u(x) = e”. 
Applying the theorems (1-7), the differential transformation of the problem ( 3. 9 
) can be determined as 


(3. 11) U(k4 


1)"U(k — ke 





13) = 


1, 


JU (ka = kı) 





nms ! 


kı! 


Using Eq. ( 2. 1 ), the boundary conditions ( 3. 10 ) at zo = 0, can be transformed 


as 


(3. 12) 


(3. 13) 


II 
i 
G 
~ 
jæi 
w 
II 
q 
e 
— 
N 
atar 
Sj= 


2-0 U (k) = 


2)(k — 3)U(k) = e, 





2)(k — 3)(k — 4)U (k) = 


€, 
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where n is a sufficiently large integer. Using the recurrence relation ( 3. 11 ), the 
transformed boundary conditions ( 3. 12 ) and the inverse differential transforma- 
tion ( 2. 3 ), the following series solution upto O(a'*) can be determined 








r2 r3 xt gp? zê 
= 1 ! ! ! Ag’ + Bz? 
u(x) +2 25 24 + 190 799 + e+ Ba? + Cx 
13 
14 +Da! + Bart! + Fr? 4 a r 
id) = "+ Sea7020800 | OF): 


where, according to the definition ( 2. 1 ) 


uD ul) u) 
A220 E SOW 78) C= yg), 
(12). 





(3. 15) yo) PAGED) oF a2) 

p=+,© =U(10), E = Ke =v), F=4,M =U 
From the transformed boundary conditions ( 3. 13 ), the following system of linear 
equations in terms of A, B and C can be determined 


16925388481 | A+ B+C+D+E+F=e, 


11! 





1301287681 = 
“arapoieou. T 14 t 8B +9C +10D +11E +12F =e, 


108108001 y 
30916800 + 424 + 56B + 72C + 90D + 110E + 132F =e, 


(3. 16) 


7 | | | — 
oer eens + 210A + 336B + 504C1 + 720D1 + 990F + 1320G = e, 








seta + 840A + 1680B + 3024C1 + 5040D1 + 7920F + 11880G = e, 


SiL + 2520A + 6720B + 15120C1 + 30240D1 + 55440F + 95040G = e. 
The solution of the above system, yields 
A = 0.00019841261021669805, B = 0.000024802098855594142, 
C = 2.75451301223356 x 107, D = 2.7708206549274656 x 1077, 
E = 2.4060411409761992 x 1078, F = 2.3783375635762833 x 107°. 
The series solution can, thus, be written as 








x? z3 x4 x zê 
= 1x4 + 0.000198413a" + 0. 24802128 
u(x) l+e BE ae CTRET 0.000198413a" + 0.00002480212 
+2.75451 x 1076x? + 2.77082 x 1077! + 2.40604 x 1078x! 
13 
zi 42.37834 x 107941? + ——* _ 4. Q("4). 
(3:17) 37834 x 107°?x a e o (x14) 





The approximate solutions of Example 3.2 obtained with the help of DTM are 
compared with the results of the VIM [1] in Table 2. From the numerical results, 
it is clear that the DTM is more efficient and accurate. By increasing the order 
of approximation more accuracy can be obtained. The graphical comparison of 
exact and approximate solutions is shown in Figure 3. The solid line represents the 
curve corresponding to the exact solution whereas the dotted line corresponds to 
the approximate solution. In Figure 4, absolute errors are plotted. 
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FIGURE 3. Comparison of 


the approximate solution 
with the exact solution for FIGURE 4. Absolute errors 
problem 3.2. Dotted line: for problem 3.2. 


approximate solution, solid 


line: the exact solution. 


Conclusion In this paper, the differential transformation method has been 
applied to obtain the numerical solution of linear and nonlinear thirteenth order 
boundary value problems. The present method has been applied in a direct way 
without using linearization, discretization, or perturbation. By increasing the order 
of approximation more accuracy can be obtained. Comparison of the numerical 
results with the existing technique [1] shows that the present method is more 


accurate. 


TABLE 1. Comparison of numerical results for Example 3.1 


















































x Exact Approximate | Absolute Error | Absolute Error 
solution | series solution | present method [1] 

0.0 | 1.0000 1.0000 0.0000 0.0000 

0.1 | 1.09484 1.09484 2.22045E-16 3.88578E-15 
0.2 | 1.17874 1.17874 0.0000 1.46216E-13 
0.3 | 1.25086 1.25086 2.22045E-15 8.80518E-13 
0.4 | 1.31048 1.31048 6.66134E-15 2.35822E-12 
0.5 | 1.35701 1.35701 1.11022E-14 3.8014E-12 
0.6 | 1.38998 1.38998 1.04361E-14 5.14766E-11 
0.7 | 1.40906 1.40906 5.32907E-15 1.56224E-11 
0.8 | 1.41406 1.41406 8.88178E-16 8.99409E-11 
0.9 | 1.40494 1.40494 0.0000 4.70031E-10 
1.0 | 1.38177 1.38177 0.0000 2.06386E-09 
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TABLE 2. Comparison of numerical results for Example 3.2 
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x Exact Approximate | Absolute Error | Absolute Error 
solution | series solution | present method [1] 
0.0 | 1.0000 1.0000 0.0000 0.0000 
0.1 | 1.10517 1.10517 4.44089E-16 4.17444F-14 
0.2 | 1.2214 1.2214 4.44089E-16 2.64144E-12 
0.3 | 1.34986 1.34986 2.44249E-15 2.99314E-11 
0.4 | 1.49182 1.49182 7.32747E-15 1.67101E-10 
0.5 | 1.64872 1.64872 1.22125E-14 6.30955E-10 
0.6 | 1.82212 1.82212 1.11022E-14 1.84757E-09 
0.7 | 2.01375 2.01375 5.77316E-15 4.47866E-09 
0.8 | 2.22554 2.22554 1.77636E-15 9.21592E-09 
0.9 | 2.4596 2.4596 8.88178E-16 1.58906E-08 
1.0 | 2.71828 2.71828 0.0000 2.09057E-08 
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